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Abstract 

I present the exact energy eigenstates and eigenvalues of a quantum many-body system of bosons 
on non-commutative space and in a harmonic oszillator confining potential at the selfdual point. I also 
argue that this exactly solvable system is a prototype model which provides a generalization of mean 
field theory taking into account non-trivial correlations which are peculiar to boson systems in two 
space dimensions and relevant in condensed matter physics. The prologue and epilogue contain a few 
' remarks to relate my main story to recent developments in non-commutative quantum field theory and 

an addendum to our previous work together with Szabo and Zarembo on this latter subject. 

O 

(N 

^ , Prologue 

' Recently the renormalizahility of certain non- commutative (NC) quantum field theory (QFT) models was 

established in important work by Grosse and Wulkenhaar and the Orsay group [1~5]; see Ref. [6] for 
fT^ ' further references and a lucid discussion of the significance of these results. It seems that all examples of 

renormalizeable such models share the duality property that their actions have the same form in Fourier- 
' (—( '' as in position space [7]. In this contribution I will discuss the following models possessing this latter 

O |. property [8-11]: 

i ■ n 

: s 



\ = [ d^"x {<i>''{x)[a{-id- B ■ x)^ +a{-id + B ■ xY - H]<i>{x) +g^^ i.^i.'i)^ *'i>{x)) . (1) 



[Notation: I denote points in M^", 2n = 2, 4, . . ., by x = {x^ , ■ ■ ■ , a;^"); ^'^^\x) represents a boson field to 
be specified in more detail below; I write 



2 2n 



0^: (-ia±B-a;)2 = ^ f-i— ±(B.x% 1 with (B.x)^ = }_^B^,x^ (2) 

oo 

I and B = some fixed skewsymmetric and invertible 2?! x 2?! matrix; a, a > 0, fi and g are real 

i parameters such that cr -f ct > 0; ★ is the well-known Groenewold-Moyal product (see e.g. [12] for review) 

characterized by another skewsymmetric 2n x 2n matrix = {& fi:j)f[^^=i as follows, 

O ■ x'^^x'' -x"ir x" = -216'^''; (3) 

_ the dagger indicates complex conjugation and the Hilbert space adjoint; N are the positive integers.] / 

■ denote the same mathematical expression by two different symbols since it can be either interpreted as the 

H I Hamiltonian (Ti) of a quantum many-body system on 2n dimensional space [8, 9] or, alternatively, as an 

action (S) of a NC QFT model on 2n dimensional Euclidean spacetime [7, 10, 11]. My discussion will 
be mainly restricted to the special case where BO — I (i.e. is the inverse of the matrix B) where these 
models are exactly solvable [8, 11]. One of my aims is to present results which have remained unpublished 
up to now, another to add a few remarks to our previous publications on this exactly solvable NC QFT 
model together with Szabo and Zarembo [10, 11] to point out an interesting alternative interpretation of our 
results, and to emphasis that an interesting problem concerning this model still remains to be solved. The 
new results are on the exactly solvable quantum many-body system of bosons and presented in section.s[^^ 
and my remarks on the NC QFT model are contained in the epilogue at the end. 

Before going into my main story I shortly recall the simplification arising at BO = / [8] (in this discussion 
I will refer to the mathematical object in {ip as Hamiltonian TL, but everything I say applies word-by- 
word also to its interpretation as action S). Obviously, the Hamiltonian in {Ip is the sum of two term, 
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7i = Ti.Q + Ti.int, where Tio and Timt are the quadratic and quartic parts in the fields respectively. As is 

well-known, in many standard field theory models one can expand the fields in a basis such that either the 
quadratic or the quartic part of the Hamiltonian becomes simple, but in general it is not possible to make 
both parts simple in the same basis, and this is one main reason why, in general, field theory models are 
computationally challenging (typically Ti.Q is simple in Fourier space and Tiint in position space). However, 
for the model in at the special point BO — I there exists a basis in which both, TLq and TLmt, oltc simple. 
To be specific: this latter basis is given by the common eigenf unctions of the differential operators 

in 0j labeled by two sets of integer vectors £,m € N", and by expanding the fields in this basis 

$(^)-E^^»'^L(^)' *^(^) = E^L'^^™(^) (4) 

£,rn 

the Hamiltonian acquires the following remarkably simple form [8, 11], 



n 

s 



The parameters Eg and E^-^ are proportional to the eigenvalues of the operators in 0j and given by 

n ^ n ^ 

E,^Y.^a\B,\{l,~-), Er^=Y.^d\B,\{m,~-) (6) 



2' 

j=i j=i 

with £ = {£i, £2, . . . , £n) cind similarly for m; \Bj \ are eigenvalues of the matrix \J B^ (see [11] for a precise 
statement), and 

9= ~' . (7) 

It is interesting to note that the model in ^ can be written in the following matrix form, 

= Trace(£;AU + EAA^ - fiA^ A + g{A^Af ) (8) 



n 

s 

where A, E and E above stands for the infinite matrices with matrix elements Aim, EpS^ra, md EmS(rm 
respectively, the matrix adjungation is defined such that {A^)^^ — Aj^^, and matrix multiplication is un- 
derstood. 

1 Introduction 

Interacting boson systems have been of interest in theoretical physics since the early days of quantum 
physics, and a recent increased interest in this subject was triggered by remarkable experimental progress 
to realize and study the Bose-Einstein condensation; see e.g. [13] for a recent text book in this topic. I 
believe that these developments provide a good additional motivation for studying the NC quantum many- 
body Hamiltonian Ti. in as I will argue in more detail in my first remark in section [H this Hamiltonian 
defines a prototype model which allows to study a particular type of correlations and its effect on the 
Bose-Einstein condensation in an exact solution. 

In the main part of this paper I thus interpret the model in ([T]) as Hamiltonian Ti. of bosons moving 
on 2n dimensional space M^" and interacting with a particular four point interaction. I will show that 
this model is exactly solvable in the sense that all its energy eigenstates and eigenvalues can be computed 
explicitly. As will be seen, this exact solution provides an example of a correlated boson system. To simplify 
notation and to allow for a simple physical interpretation I restrict my discussion to the case 2n = 2 and 
(7 = (T = 1/4, but my results can be straightforwardly generalized to 2ri > 2 and general parameter values. 
The parameter fi corresponds to the chemical potential, and, for my purposes, one can assign to it any 
convenient value. 

I mention in passing that the fermion variant of this model was introduced and analyzed in [8,9] but, 
to my knowledge, the boson story presented here has not appeared in the literature before. 

The plan of the rest of this paper is as follows. In section [2] I give a precise definition of the quantum 
many-body model, and in section [3] I present its solution. Section 2] contains various remarks, and, in 
particular, I explain there why I believe that this model is a prototype model for interacting bosons and 
relevant in condensed matter physics. 
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2 Definition of the model 



I consider the quantum many-body system defined by the Hamiltonian in llj where the boson fields <i>^^-' (x) 
are operators acting on a boson Fock space T defined by the usual canonical commutator relations and a 
normalized vacuum state annihilated by all operators ^{x)\ see e.g. [14]. Expanding the fields as in ^ 
these latter relations are equivalent to 

for all £, TO, £' , m', as usual. Choosing B9 ~ I this Hamiltonian can be written as 

i,m £,m.i' ,m' 

where I found it convenient to normal order the interaction term Tiint (this corresponds to a renormalization 
of /i which can be ignored) and rename ~ ji to Eg. To simplify my discussion I set 2n = 2, a = a = 1/4 
and /i = \B\ so that ^. m g N and 

E, + E„,^\B\{l + m~2). (11) 

Then the model describes interacting bosons on two dimensional space and confined by harmonic oszillator 
potential, with £ — 1 and to — 1 the usual harmonic oscillator quantum numbers. A useful alternative 
interpretation of the quantum numbers I and m is as x- and y components of (quasi-)momcnta of bosons in 
two dimensions. In this latter interpretation one is interested in other dispersion relations like Ei + E,-a cx 
(^■^ + TO^) and m nmning also over negative integers (see the first remark in section|4]), and it is therefore 
important to note that such changes do not affect the exact solubility of the model. Then the interaction 
term describes two body scattering processes where two bosons with initial momenta (£,to) and {£',m') 
exchange the y-components of their momenta while the x-components remain the same, or vice versa [9]. 



3 Exact solution 

I now discuss how to construct exact energy eigenstates of this model. For that I consider the quasi-free 
states (the normalization of the eigenstates will be ignored) 

for fixed quantum numbers £j and rrij in N, with N an arbitrary fixed non-negative integer. One can 
interpret this as a state containing N bosons with momenta {£j , rrij ) . Each such state is an eigenstate of 
the quadratic part Tip of the Hamiltonian, and the corresponding eigenvalue is 

N 

So = Y.{E,^+E^,)- (13) 
i=i 

I will refer to this as kinetic energy. 

It is important to note that the permutation group Sn of N elements has a natural action on these 
states T] as follows, 

for all P e Sm, and that all these states Prj are degenerate eigenstates of Tio- Moreover, one can show that 
the action of the interaction part of the Hamiltonian TLint on such a state rj is 

TiintV = 2g ^ TjkT] (15) 
i<j<fc<w 

where Tjk € Sn is the transposition which interchanges j and k and leaves all other integers 1, 2, . . . , iV the 
same. One can interpret Tjk as the operator exchanging the y-components of the momenta of the j-th and 
the fc-th boson leaving the x-components the same, or vice versa. Obviously this implies that all eigenstates 
of TL are of the form 

= ^ apPr/ (16) 
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for some 77 and certain coefficients ap to be determined, f will get back to the problem of how to construct 
all these eigenstates and corresponding eigenvalues further below. 

For now I consider particular such eigenstates which can be obtained by elementary methods and which 
include the groundstates in the weak- and strong coupling limits. These eigenstates are given by 

where (+)^ is always 1 and (— )^ = f for even and —1 for odd permutations P, respectively. To see that 
these are eigenstates we note that Tjk'^± = ±5*^, which implies is an exact eigenstate of Tiint with 
eigenvalue zLgN[N — 1) (since X/i<j<fe<7v ^ 

7V(iV- l)/2), and thus 
Hv^ = {£o±gN{N-l))7j^. (18) 

It is interesting to note that the state rj^ has a fermion-like character and, as discussed below, this implies 
a strong variant of the Pauli exclusion principle which will play an important role. As will be shown further 
below, the states r]^ are extremal in the sense that they have the largest possible interaction energies. 

In particular, for g < the groundstate of the model at fixed particle number N is the state 77+ such 
that the kinetic energy £0 assumes its smallest possible value. It is easy to see that the state in (fT^ with 
the minimum kinetic energy is 

m = {Al.fn, (19) 

i.e. all bosons are in the same one-particle state with momentum {£,m) — (1,1). Note that rji is the 
well-known Bosc-Einstein condensate (BEC) groundstate of the non-interacting system {g = 0). In fact, 
this state is the groundstate for all g < Q (this is true since rj^ equals 771 up to a constant), and it is easy 
to see that the corresponding groundstate energy is 

£i=gN{N-l). (20) 

For g > the states rj^ have a large interaction energy, and for sufficiently large g > the groundstate 
of the model should be the state ri~ with 77 such that the kinetic energy Sq is minimal. It is important to 
note that one now cannot take as r] the state 7/1 in (|19p since 77^" vanishes. More generally, the following 
strong variant of the Pauli exclusion principle holds true: The state 77^ in Jj?] ) is non-zero only if all the 
X- and all the y-components £j and rrij of the boson momenta in the state rj in hlS]) are distinct^ [Proof: 
Consider a state 77 in (jl2p such that Ij = £k and/or nij = for some j < k. This implies Tjkr] — rj, 
but then rj^ = ^p{ )^Pri = ^p{ )^ PTjkrj = —rj^, and thus rj^ = 0.] A state whose momenta are all 
distinct and which has the lowest possible kinetic energy is 

^72 =4,i4,2---^]v,jv^' (21) 

and thus 

is the groundstate of the model in the strong coupling limit. The corresponding groundstate energy is 

£2^\B\N{N -1) ^ gN{N -1) (23) 

(since J^Zii^ - 1) = ^(^ - l)/2)- As discussed below, 772 is actually the groundstate of the model not 
only in the strong coupling limit but for all g > 

I now discuss the problem of finding the groundstate for intermediate coupling values. Note that, for 
fixed ?7 in (fT^ . the states Prj in (fn|. P g Sn, span a subspace !Frj of the boson Fock space J^. The 
dimension of this subspace is < A^!, and it is A^! if and only if all the x- and y-components £j and mj of 
the bosons in the state rj are distinct. It is important to note that defines a representation of the 
permutation group Sn on J^jp and this representation is, in general, reducible. Moreover, the operator 

Cat = ^ Tjk 

l<j<k<N 

^The standard Pauli principles for formions is weaker since it only requires that all the momenta (£j,mj) are distinct. 
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appearing in p5p commutes with aU permutations P G Sn, and it is therefore a constant in each irre- 
ducible representation (irrep) of Sn- Since on J-^ the kinetic energy Tio is constant and the interaction 
Ti-int proportional to Cn, the problem of constructing eigenstates of TL is equivalent to decomposing the 
representation of Sn on J-^i described above into irreps. This is a classical problem solved in group theory; 
see e.g. Chapter IV in [15]: the irreps of Sn can be labeled by partitions A of N, i.e. A = (Ai, A2, . . . , A^-) 
with integers Xj such that 

K 

Ai > A2 > • • • > Aa' > 0, = ^' (24) 

and the value of Cat in an irrep A is Cat = J2f=i[^^ji^j + 1) ^ jAj]; see e.g. Eq. (4-3) in [16]. Moreover, the 
states in an irrep A can be obtained by applying to states 77 the so-called Young symmetrizer [15] denoted 
by y'^. One thus concludes that the energy eigenstates of this model are F'^ry with the corresponding 
eigenvalues 

N K 

^ = E (^^. + ) + E 9>^^ (Aj + 1 - 2j) . (25) 

The states 77+ and 77^ in (fT7|) correspond to the special cases A = {N) and A = (1, 1, . . . , 1) = (1^), respec- 
tively. In principle this gives all eigenstates and eigenvalues of the model. There is, however, an important 
complication: as seen in the previous section for the special case A = (1^), if there are degeneracies many 
of the eigenstates Y^'^ri vanish. To find the groundstate of the model we therefore must determine the state 
r] in p2p of minimal kinetic energy and such that, for a fixed partition A, Y^rj is non-zero. This problem 
has the following solution, 

r;=(4i)^n4,2)'^---(^U)'"^^, (26) 
and the smallest possible energy eigenvalue in an irreps A is therefore 

K 

= E (^^ + + 9(^1 + 1 - 2^)) • (27) 

i=i 

One can determine the groundstate of the model by finding the partition A of iV which minimizes the 
energy in ([?7|) . The solution of this problem depends on the dispersion relation Ei + E„i. Using the one 
in (fTT|) one finds A = {N) for g < and A = (1^) for g > \B\, which confirms that the states in (fT^ and 
([21]) are the groundstates for g < and g > \B\, respectively. In the intermediate regime < g < \B\ the 
groundstate is given by a partition approximated by 

/ O AT I U\ 

(28) 

where "~" means that the l.h.s. is the non-negative integer closest to the r.h.s., and this approximation 
becomes exact in the limit when the boson number N becomes infinite. 

I finally note that one can prove that the eigenstates of this model are, in general, correlated by finding 
one non-zero connected 4-point correlation function. As an example I consider the normalized strong 
coupling groundstate for iV = 2: 

* = i= (4,14,2 -4,24i)"' (29) 

which supports the following non-trivial connected 4-point correlation function, 

(*, 4,24,1^1,1^2,2*) - (*,4,2A,1*)(*4,1^2,2*) - (*,4,2^,2*)(*4aA,l*) = (30) 




with (•, •) the inner product in the boson Fock space. 
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4 Concluding remarks 



A key problem in theoretical physics is to do reliable computations in quantum models with interac- 
tions so large that perturbation theory does not apply. One well-known and often successful strategy 
in this context is mean field theory. It is interesting to note that one approach to mean field theory is 
to truncate the interaction in the model under consideration and only keep the so-called Hartree- and 
Fock terms, which typically leads to an exactly soluble model whose solution is equivalent to mean 
field theory of the original model; see [9] for a discussion of this in the context of interacting fermion 
systems. Mean field theory does not take into account correlations, and the latter are believed to be 
particularly important in two spatial dimensions (2D). I propose that the Hamiltonian in I11U\} defines 
a prototype model allowing to study important 2D correlations in an exact solution. To motivate this 
I consider the following standard 2D boson model 

^ = E^^^('')^W + Z^ E ^k,+k„k3+k,6nki)&t(k2)6(k3)6(k4) (31) 

k ki,k2,k3,k4 

with the boson mass M > and coupling parameter U > 0. The boson operators fe^^'(k) are labeled 
by 2D momenta 

27r TT 
k= {kx,ky), kx,ky G -pZ such that \k.j;,y\ < — (32) 

and obey the usual relations, [6(k), 5t(k')] ~ (5k, k' etc. The parameters L ^ a > correspond to 
the system size (L) and a lattice constant (a) and provide a IR and UV cutoff for the model. The 
interaction term in this Hamiltonian comes from a local interaction in position space and describes 
scattering processes where two bosons with momenta ks and k4 are scattered into states with momenta 
ki and k2, and the model is complicated since all possible such scattering processes occur with equal 
strength and are restricted only by overall momentum conservation. The Hartree- and Fock terms 
correspond to the scattering terms where ki = k4, k2 — ka and ki = ka, k2 = ka, and they are 
(essentially) trivial for this model in the sense that they only add an energy ~ AgN^ and do not 
(much) affect the groundstate. Note that the interaction contains also the scattering terms where 

(ki)x = (^4)2,, {k2)x = ik3)x, {ki)y = {h)y, (^2)?, = (^4)?, (33) 

and similar terms with x and y interchanged. These scattering terms are Hartree-like in the x- 
and Fock- like in the y-component of the momenta and vice versa, and they are peculiar to 2D. // 
one restricts the interaction terms in the Hamiltonian in i31\) and only includes these latter mixed 
Hartree-Fock terms one obtains exactly a Hamiltonian as in ilOp with 

E, + E,^ = ^(^^^\f + m% (34) 

and integers i,m such that \i\, \m\ < L/(2a). As mentioned, this latter truncation is very similar to 
a successful method to derive useful mean field theories for interacting fermion systems, and I thus 
regard the model in (fTU|) as a generalized mean field model. The exact solution of this model above 
does not rely on the form of Ei + Em (except for the groundstate, of course), and, as argued below, 
this model describes interesting "physics" which cannot be accounted for in mean field theory. 

It is interesting to note how the character of the groundstate of the model changes with increasing 
coupling constant g: for .g = one has the standard BEG groundstate in p9|l where all bosons are in 
the same one particle state {£j,mj) = (1,1). As the coupling increases it becomes more favorable to 
reduce the degeneracies and thus the number of bosons in the BEG, and one finds a distribution of 
the momenta as described by the partition in and a correlated groundstate. Moreover, the BEG 
condensate in the ground state for 2|i?|/(A^ -f 2) < 5 < is 



(4,1 ^ ^ ^2iV (^M „ 1^ , (35) 

and it becomes 1 for g > \B\ where the groundstate becomes maximally correlated. Moreover, 
as demonstrated in pop above in a simple example, one can construct non-trivial 4-point Green's 
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functions for the model to prove that its groundstate is, in general, correlated and thus not accessible 
by mean field theory. I hope that these remarks are sufficient to convince the reader that the "physics" 
of this model is non-trivial and interesting. I plan to present a more detailed discussion elsewhere. 

3. It is important to note that the model in (jlOp describes a stable system only in the parameter 
regime < g < \B\ since otherwise the groundstate energy can be decreased by increasing the 
particle number N to infinity (this follows from ((20)) and ([25)) ). In my interpretation of this model 
as generalized Hartree-Fock model the instability for g > \B\ is removed by the Hartree-Fock energy 
~ AgN^ which should also be included. 

4. Obviously nearly everything I wrote in the previous section can be immediately generalized to 2n > 2 
and other values for a and ct, and the only change will be the solution to the problem to minimize 
the energy in ([77)) . 

5. In this paper I only computed the groundstate of the model and demonstrated how to compute the 
other energy eigenstates and eigenvalues. Obviously it would be interesting to compute also other 
quantities, like Green's functions and the partition function. 



Epilog 

The two models in {I)) are closely related: the NC QFT model defined by S can be obtained as infinite 
temperature limit of the quantum many-body model TL. Indeed, one can write the generating function for 
the Green's functions of the latter model as matrix path integral (where the integration variables A^^i^iT) 
are periodic functions of the Matsubara time r £ [0,/3] with [3 the inverse temperature; see e.g. [14-]), ^.i^d 
the functional integral defining the NC QFT model [11] can be obtained as a limit /3 — > from that. Thus 
the model Ti in (QP defines a 2n + 1 dimensional QFT. It would be interesting to use this relation to defer 
from my results on the latter model results for the former model. 

The NC QFT model S in (OJ) has interesting and non-trivial QFT divergences which one has to treat 
by regularization and renormalization. I believe that this models provides an interesting example where the 
role of such divergences can be studied in detail and by exact and explicit results beyond perturbation theory, 
and that our previous results on this [10, 11] are only a first step in this direction: as I will argue below, 
there are other QFT limits than the ones studied in this latter work, and one of these other limits is more 
interesting and more difficult than the others. The explicit solution of the model in this latter limit is a 
challenging but doable project for the future. 

To be more specific: The natural regularization for the NC QFT model S in is to restrict the fields 
A'iH to {£i,£2,---,in) such that 

gj = l,2,...,L<oo (36) 

and similarly for m. Then the fields Agm can be naturally interpreted as components of a^ N x N matrix A 
with N = L". With that the functional integral defining the NC QFT model becomes a well-defined integral 
over , and the non-trivial task is to find a dependence of the model parameters a, a, —ji, and g on the 
cut-off parameter N such that the limit N ~^ oo is well-defined an non-trivial. 

In [11] we studied two such limits for the case a = which we called IR- and UV limit: The IR limit 
corresponds to the following scaling of parameters, 

-, gren j-, Bren /nr,\ 

^ = 1, a=—, B=j^, fi^tiren (37) 

where the parameters with the subscript "ren" (short for "renormalized" ) are independent of N . The results 
for the Creen's functions in this limit can be found in [11]. I only mentioned here that B for TV — > oo 
leads to a 2-point Creen's function which is translational invariant, G{x;y) — G{x — y;0), and all the higher 
Creen's functions are trivial. In the IR limit the duality symmetry of this model [7] is broken, which implies 
the existence of a dual limit where B ^ oo as N —^ oo and with Creen's function obtained from the ones in 
the IR limit by a duality transformation. The 2-point Creen's function in this latter UV limit is non-trivial 
and ultra-local, G{x;y) oc S'^"{x — y), and the higher Creen's functions are again trivial [11]. 



'Note that the symbol A'^ in the following and in sections 11141 have different meanings! 
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Now comes my addendum to [11]: It is possible to get a third limit in which the above mentioned 
duality symmetry is not broken as follows: rather then keeping a constant and scaling B like A^^^/" one 
can scale a like N~^/^ and keep B constant in the limit N — > oo: 

Cren Qren j-, j-, /oo^ 

The non-trivial scaling of ct can be interpreted as multiplicative regularization. It is easy to deduce from 
the results in [11] the 2-point Green's function in this third limit, 

G{x;y) = Y.(<nAJ^imiMmiy) (39) 

i,m 

with (^|„j^£„j) depending only on £ and computed explicitly in [11], section 4.2; the higher Green's functions 
are again trivial. 

I finally would like to emphasis that the limits described above are restricted to the case a — 0, and 
for a > there should exist another limit leading to non-trivial higher Green's function and which should 
describe a non-trivial fixed point of the renormalization group of the models in (OJ). / expect that this latter 
limit is the one studied in a closely related model in [17, 18] (this latter model is similar to ours for a = a). 
It is certainly not easy but, as I believe, possible and very desirable to compute the Green's functions of the 
model in this QFT limit explicitly. In this context I should mention the non-perturbative renormalization 
of NG (l>^-theory which was recently established by Grosse and Steinacker [19]. 

Note added: The model ?i in ([T]) for 2n = 2, cr = ct = 1/2, and at finite temperature was recently studied 
in [20]. 
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